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SUMMARY 


A matrix solution Tor the spanwise bending moments and deflections 
of a torsionally rigid rotor blade subjected to an arbitrary loading is 
presented. The method includes the cantilever, teetering, and hinged 
blades in hovering and in steady forward flight. The method is compara- 
tively short, involves only standard matrix procedures, and does not 
require that the mode shapes or natural frequencies be known. 


INTRODUCTION 


Numerous methods are available for calculating the uncoupled span- 
wise bending moments and deflections of rotor blades in steady unaccel- 
erated flight. Some of these methods (refs. 1 and 2) involve tabular 
solutions of the differential equation of blade deformation. These 
tabular procedures are very lengthy, particularly when numerous loading 
conditions and rotor speeds are being investigated. Other methods 
(refs. 3 and 4) require the calculation of the natural-mode shapes and 
frequencies of the blade. The calculations of mode shape and frequency 
are laborious and either must be repeated for each rotor speed or an 
approximate correction must be applied to account for the effects of 
rotor speed. A matrix method which avoids some of the difficulties of 
the tabular and modal solutions is presented in reference 5; however, 
tabular solutions of the equation for blade bending moment under quasi- 
static conditions are required before the method can be applied. 

In the present method the differential equation for blade bending 
moment is solved entirely by matrix procedures. The method permits the 
determination of the bending moments and deflections directly without 
preliminary quasi-static, mode-shape, or natural-frequency calculations 
and, as a result, the method is shorter and much work is eliminated. 

The method is adaptable to cantilever, semirigid (teetering), and fully 
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articulated (hinged) "blades in hovering and in steady forward flight and 
is applicable to design calculations and experimental data analyses. 

The differential equation for the blade bending moment is composed 
of matrix expressions for the centrifugal inertia, vertical inertia, and 
damping loadings, and the contribution of each term to the blade loads 
and moments may be determined. The effect of blade structural damping, 
which has been omitted in other methods, is included in the present 
method. Comparisons are presented between the bending moments obtained 
by this method and the method of reference 3 for the first harmonic 
loading on a fully articulated blade and the second harmonic loading on 
a cantilever blade. 


SYMBOLS 


Some of the symbols listed in this section are illustrated in 
figure 1. 

A blade -element average aerodynamic damping coefficient, 

- ± pac Ar r, ll3 - se - c2 
2 ' In. 

a blade -element lift -curve slope, per radian 

c blade chord (at center of blade element) , in. 

El blade bending stiffness, lb-in.^ 

F centrifugal force acting on mass of blade element divided 

by lb-sec^ 

g structural damping coefficient 

/ 

h perpendicular distance from axis of rotation to flapping 

hinge, in. 

i imaginary component of complex number, /-I 

l blade running load, lb/in. 

L load on blade element, lb 

bending moment, lb -in. 


M 
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Subscripts: 

a 

C 

d 

H 

i 


mass of blade element, — 

in. 

number of blade elements for extended matrices 
distance from center of hub to blade tip, in. 

spanwise distance from center of hub (measured along 
undeflected blade) , in. 

width of blade element, in. 

forward velocity, in. /sec 

elements of an influence coefficient matrix (see table III) 


blade deflection, measured from a plane of rotation per- 
pendicular to shaft axis, in. 

rigid-blade flapping angle, measured from a plane of 
rotation perpendicular to shaft axis, radians 

angle of flexible blade at hinge, measured from a plane 
of rotation perpendicular to shaft axis, radians 

blade slope due to blade bending, dz/dr, radians 


rotor angular velocity, radians/sec 


angular frequency of applied loading, radians/sec 


air density. 


lb-sec* 1 


xn. 


blade azimuth angle (measured in direction of blade rota- 
tion from downwind position) , radians 


aerodynamic 
cosine component 
damping 

due to forces parallel to plane of rotation 
inertia 
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L pertaining to a load station 

r pertaining to a spanwise station 

s structural 

S sine component 

1,2, . . .n station number or number of harmonic 

V due to forces perpendicular to the plane of rotation 

w blade weight 

z deflection 

f centrifugal loading 

A dot over a symbol indicates first derivative with respect to time. 
Two dots indicate second derivative with respect to time. 


Matrix notation: 

t ] 


[ i 


rectangular or square matrix 
diagonal matrix 


{ 

[ 

L 


} 


l 4 

J 


column matrix 

inverse of a square matrix 
row matrix 


METHOD 


The differential equation for the structural bending moments of 
hinged, teetering, and cantilever rotor blades is set up in matrix form. 
Matrix expressions for the blade moment due to the various blade loadings 
are derived and added to form the blade -moment equation. The differen- 
tial equation is first derived in a generalized form and is later adapted 
to (l) the case where the rigid-blade aerodynamic loads are known from 
design calculations and the flexible -blade structural moments or deflec- 
tions, or both, are desired; ( 2 ) the case where the blade structural 
moments are known from strain -gage tests and the total aerodynamic 
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moment, including aerodynamic damping, is desired; and (3) the case 
where the aerodynamic loads have heen measured and the structural moments 
are desired. The procedure for obtaining the rotating-blade character- 
istic matrix is also presented so that the natural modes and frequencies 
may be determined. 


Blade Loadings 

In the calculation of spanwise bending moments, a torsionally rigid 
rotor blade or its equivalent (an elastic blade for which the section 
moments are zero and the aerodynamic center, elastic axis, and center of 
gravity coincide) may be treated as a simple beam subject to numerous 
superimposed loadings. If chordwise effects are excluded, the blade in 
the forward flight condition is subjected to the following distributed 
loads : 


(1) The aerodynamic load L a (excluding the aerodynamic damping 
effect which is directly proportional to blade flexure plus flapping 
deflection velocity; see item 5 ) 

( 2 ) The centrifugal load F 
( 5 ) The blade weight 

(4) A normal inertia loading due to vertical accelerations Lj_y 

( 5 ) An aerodynamic damping load (defined by eq. ( 5 )), directly 

proportional to blade deflection velocity due to both flapping and 
flexure 

( 6 ) The blade structural damping load L^ s 

The combined moments of these six loadings are balanced by the blade 
internal moment due to blade bending resistance. 


Derivation of Blade Moment Equations 

The blade geometry, the forces acting on an infinitesimal ele- 
ment of the blade under forward-flight conditions, and the resulting 
force system are shown in figure 1. From figure 1 it is seen that the 
bending moment M rn at any station r n due to the loads acting on the 
blade is given by: 
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nR o R-' 

M rn = / M r “ r n)^ r - / *h( z - *n)* r C 1 ) 

J r n J r n 

The matrix equivalent of equation (l) may be written in the form 
{M} = [r]{L v j - n 2 [F]{z} (2) 

Equation (2) corresponds to a step-diagram representation of the 
Zy and Zg loadings of equation (l). The Qr] matrix is an integrating 

matrix with values that depend on the stations and integration method 
used* The QfJ matrix incorporates the values of the horizontal loads 
{ l 1h} (dimensional with respect to rotor speed) as well as integrating 

numbers. The £r] and [Fj matrices are derived in the appendix. 

Substituting the five components of |Lyj* enumerated previously 

into equation (2) yields the following bending -moment equation of the 
combined loadings: 

{Mj.} = H{L a | + [r]^} + [r]{L da J + [r] {L iV } + [r] {L ds } - fi 2 [F]{z} 

( 3 ) 

Equation (3) is the general form of the equation for blade bending 
moment used in this study and pertains to hinged and cantilever blades 
in hovering and forward flight. In order to obtain the blade bending 
moment for any specific loading condition, the matrix expressions for 
the various loads acting are substituted into equation (3). 

Loading matrices . - The expressions for the various component loads 
in equation (3) become quite simple if equation (3) is restricted to the 
expression of the moments due to loadings of a single frequency. The 
total aerodynamic loads on the blade may be broken down into components 
which are of a sinusoidal nature, and the standard procedure is to 
develop equation (3) for the moment due to a single frequency. Super- 
position may then be used to obtain the total moments due to all of the 
harmonic components on the blade. The loading matrices which appear in 
equation (3) are given by the three following equations. 

For a particular frequency ao of the applied sinusoidal load, the 
vertical inertia loads |L^y| are given by 



NACA TN 4304 


7 


{l 1v } = co 2 ] (4) 

where the elements of the £m ^ matrix are the masses of the blade 
elements . 


The aerodynamic damping loads |L^ a | are given hy 

{ L da} = ia<u M W 


where the diagonal matrix [\A^ is equal to 


- i pacr Ar 


(5) 


The elements of the £aj matrix can he deduced from steady aero- 
dynamics for the lift on a section of chord c and length Ar rotating 
at unit angular velocity and with a unit normal velocity. 


The structural damping moment 


MM 


is given hy 


Hfcds} = -ig D-.o]{m s } 


(6) 


Moments and deflections at a single frequency .- The hlade loadings 
at a particular frequency to (eqs . (4) , (5 ) > and (6)) may he substituted 
into equation ( 3 ) in order to obtain the equation for the hlade moment in 
terms of the structural moment, the aerodynamic moment, and the hlade 
deflections as 

E 1 + is] (Ms) = H { L a} + ifico[r] [a] {z} + (a^[r] [m^j {z} - ^ 2 [f] (z> 

(7) 

The deflection {z} in equation ( 7 ) is composed of the flexural 
bending and the effects of blade-root rotation. The flexural bending 
is defined in the appendix as 


M - (8) 

The total blade deflection including root rotation is, therefore, 
given by 


{z} = [Z M ] {M s } + {r - h} tan p h 


(9) 
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For small angles tan p^ « p^, an d in the equations to follow, which are 
based on assumptions of small angles, p^ is assumed equal to tan p^. 
Equation (9), therefore, becomes 



( 10 ) 


Substituting equation (10) into equation ( 7 ) yields the following 
relation between the aerodynamic loads and the structural moments for 
the flexible blade: 


(1.0 + ig) 


Mi{o} 

[0] [ 0 


« 2 [f] [[%] | {r - h}J - 0)2 [M zV ] [[Zm] j {r - h}] - 


iioG[r][A] [Z M ] | {r - h} 


|H}_ 

Ph 


- H M 


(11) 


where 

[ M zv] = MM 

The partition lines of equation .('ll) effectively separate structural 
moment effects from root slope effects, both of which act to balance the 
aerodynamic moment. If the root slope p^ is set equal to zero, the 

terms pertaining to the root slope drop out and the resulting equation 
is the equation for a cantilever -blade bending moment. The cantilever- 
blade equations are given in detail in a subsequent section. 


Equation (ll) may be inverted in order to obtain the reciprocal 
matrix relation 



(l.O + ig) 


[1.0] ko} 

T61T5- 


+ n2 w [[%]!{>• - 




• *9 
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where 

K} = M {L a } 

Equation (12) gives the structural moment and root slope p^ 

on the' rotating hinged blade in terms of the applied aerodynamic moment 
jMgJ. The equation includes the effect of the aerodynamic damping, which 
is proportional to the total deflection velocity. 

When calculations are to be made for many rotor speeds, it is advan- 
tageous to put equation (.12) in the form 



where the matrix [nJ i s evaluated for each harmonic, and 

w - [w [m \{r - »}] - f[%] [[%] - »}] - 

i §Lr][A][[z M ]j{r - h}] (14) 

The matrix [jfJ need he evaluated only for each harmonic number. 

Hovering condition .- Equation (12) , which was derived for the 
forward-flight condition, would also apply to the hovering condition if 
the terms which pertain to the blade vertical velocity and acceleration 
are omitted. Because the hovering aerodynamic loads are steady, g = 0 
and o> = 0; and the equation for hovering becomes 



Cantilever -blade equations .- If the blade-root angle is zero 

in equations (12) and (15) , the result is the equation for the cantilever 
blade in forward flights 
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N - [[1.03 + ig[l.O] + o2[f] [%] - 0>2[M zV ] [Zm] - 1^H[a][%|{m s } 

(16) 

or 

{m s } = [[l.o] + ig[l.o] + 0 2 [Fl[%] - 0)2 [M zV ] [%] - 

“MWMJ’V) (1t) 

If the calculations are to be, made for, numerous rotor speeds, equation (17) 
may be put in the same form as equation (13) . 

The cantilever -blade equation for the hovering condition is: 

{m s } = [[1.0 ] + o 2 [f] [Zm]] _1 {M a } (18) 

CALCULATION PROCEDURES 
Matrix Evaluation 

The calculation procedures presented are illustrated for the hinged 
blade. The details given pertain to the ten-segment breakdown of the 
blade load, weight, and stiffness distribution shown in figure 2. The 
procedure for either a cantilever or a teetering blade is a simplifica- 
tion of the hinged-blade case. The particular blade-section breakdown 
in the matrix tables (see figs. 2 and 3) was selected to provide more 
stations near the root of the blade where rapid changes in blade moment 
and stiffness occur. The number of blade elements used may be increased 
by the procedures given in the appendix for increasing the order of each 
matrix. However, it is believed that for most blades ten stations are 
sufficient. 

In order to apply the hinged-blade equations, it is first necessary 
to evaluate the various matrices in the equations. The hinged -blade 
equation, previously given as equation (12), is 
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'W 





(l.O + ig) 


[i.o]!{o} 


Lo J i o 


CD 


« 2 M [[Zjjj{r - h}J - 
T M zv] [M j ( r “ h }J - tea [r] [a] [[%]! {r - h}]J {Ma} 


The structural damping coefficient g is a function of the fre- 
quency, and the value used should correspond to the value of cd for 

which the equation is being evaluated. The matrix 

identity matrix where the zero, as the last element on the diagonal, 
imposes the condition of zero structural moment at the blade hinge. 

This matrix is given in table I, The matrix [f] expresses the moment 
on the blade caused by the centrifugal forces acting on the deflected 
blade. The [f] matrix is given in table II. The matrix [[Zm] | {r - h}j 

expresses the blade deflection caused by the structural moment and the 
root slope. This matrix is given in table III. The matrix [M z y| 

expresses the vertical inertia moment in terms of the blade deflection. 

The [M zV ] matrix is given in table IV. The [r] matrix is an inte- 
grating matrix which expresses the moment of a loading and is given in 
table V. The matrix [Aj expresses the aerodynamic damping load in 
terms of the blade deflection. The matrix [a] is given in table VI, 

The number of digits needed in the evaluation of the various matrices 
and in the succeeding operations for the most accurate results has not 
been rigorously established. However, limited experience has shown that 
carrying the basic masses, distances, stiffnesses, and so forth to four 
digits and, in succeeding operations, allowing the number of digits to 
accumulate to a maximum of ten gives good results » The procedure for 
filling out the given matrix tables is as follows: 

(1) Set up the identity matrix as shown in table I, which is 
already completed in its entirety and may be used as is. 

(2) Obtain the weight distribution for the blade and break it down 
into ten elements of equal length. Determine the mass of each element. 

Use the product mr for each element as described in the appendix, and 
fill out table II. 

(3) Obtain the stiffness distribution for the blade and break it 
down into the smaller elements shown in figure 2. Use these stiffness 
values and the values of Ar^ and fill out table III. 


Lo] | 0 
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(4) Use the mass values of item 2 and. fill out table IV. 

(5) Substitute the value of Ar into table V. 

(6) Use the values of the chord of each, element, the lift -curve 
slope of each element, the radius to the center of each element, and 
the values of p and Ar, and fill out table VI. 


Application to the Hinged Blade 

The equations for the hinged blade are applied to those problems 
which are believed to be the most general. The cases discussed include 
the determination of (l) the flexible -blade structural moment due to 
the rigid-blade aerodynamic loads; (2) the structural moment due to the 
measured aerodynamic loads; (3) the moment of the total aerodynamic 
loading as determined from the measured structural moment; and (4) the 
aerodynamic damping load. The procedure for obtaining the natural 
modes and frequencies of the blade is also discussed since, in some 
cases, this information may be desired. 

Structural moment due to rigid-blade aerodynamic loads .- The stand- 
ard design practice is to obtain the rigid-blade total loading from 
rigid-blade aerodynamic load and inertia balance equations . The rigid- 
blade aerodynamic and inertia loads are then applied to the flexible 
blade which is considered to flex about the rigid-blade position. This 
procedure is equivalent to the application of the rigid-blade aerodynamic 
load only (minus an aerodynamic damping component which is proportional 
to z) to the equation for the flexible blade which was derived as equa- 
tion (12) and in which the deflections are referenced to the plane of rota- 
tion perpendicular to the shaft. The equivalence of the standard design 
procedure and the procedures of the present method is proved by the 
following : 

The rigid-blade total load (aerodynamic plus inertia) which is nor- 
mally applied to the flexible blade (considered to flex about the rigid- 
blade position) is given by 


^rigid 

total 


In') - mrp - 

. /rigid 
total 


(19) 


where the terms on the right-hand side of the equation represent the 
rigid-blade aerodynamic, vertical -inertia, and centrifugal -inertia loads, 
respectively. Since the equations of the present method are written in 
terms of moment, it is useful to convert equation (19) to its moment 
equivalent. Also, it is useful to express each term in equation (19) in 
the matrix form. Thus, equation (19) becomes 
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where the terms on the right-hand side of the equation are the rigid- 
blade aerodynamic, vertical-inertia, and centrifugal-inertia moments, 
respectively. 

If the rigid-blade moments given by equation (20) are applied to the 
flexible blade, which is considered to flex about the rigid-blade posi- 
tion, the following equation results: 


Moment of rigid-blade loading 



Moment of flexible -blade loading 


(1.0 + ig){M s } + ft 2 [F] {z f lex j - “ 2 [ M zv] { z flex} " l^[r][A]/z flex | 


The terms on the right-hand side of the equation give the moments 
associated with the blade deflection -[ z pp ex j- about the rigid-blade 

position. These terms, from left to right, represent the blade -bending, 
the centrifugal -inertia, the vertical -inertia, and the aerodynamic- 
damping moment, respectively. Equation (21) implies that the blade 
moments resulting from the flexible -blade deflections must balance the 
applied rigid-blade moment. The flexible- and rigid-blade deflections 
z flex an<i z rigid 8X6 s ^ own in figure l(b) . If it is recognized that 

Zfiex Pl us z rigid e Q. ua l bo the deflection z (measured with 
respect to the plane of rotation) and if the terms of equation (21) are 
collected, then the following equation results: 

[r]{(L a ) r I = (1.0 + ig){M s } + 0 2 [F] (z> - cn 2 [M zV ] {z} - 


'rigid 

total 


i(to[r][A]|z flex } 


Adding the term -ifl^jrJfA] to both sides of equation (22), 
collecting terms, and rearranging yields 
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(1*0 + ig){M s } = [r] {L a } + iOco [r][A]{z} + <x> 2 [M zV ] {z} - 0 2 [F] {z} 

(23) 


where 

{M - (( L a ) rlgl 4 - ia»W{«w* w } W 

[ total] 


Equation ( 23 ) is recognizable as equation ( 7 ) where |L a | is 

defined as the total rigid -blade aerodynamic load minus a component 
which is proportional to the rigid-blade deflection velocity. The matrix 
[M z yj is equal to the matrix multiplication shown in 

equation ( 7 ). 


Substituting |z rigid | = {r - h}p 


in equation (24) yields 


{La} ' I( L a) rlgld l - - h}p (25) 

total J 

The flexible -blade moments, therefore, may be determined by applying 
equation ( 25 ) (the rigid-blade total aerodynamic load minus a rigid-blade 
deflection-velocity load) to the flexible -blade equation (eq. ( 12 )) which 
is referenced to the plane of rotation. 

In order to obtain the flexible -blade structural moment due to 
rigid-blade airloads: 

(1) Obtain the rigid-blade total aerodynamic load for each harmonic 
and remove the flapping effects by using equation ( 25 ). 

(2) Use the loading just obtained along with the previously deter- 
mined matrices , and obtain the moment by means of equation ( 12 ) . 

The reason the rigid-blade loads were handled in this manner was to 
keep the blade equations referenced to the plane of rotation in order to 
provide for more general applications. 

Structural moment due to measured airloads . - The procedure for the 
case of measured airloads is as follows: Since measured loads include 

all aerodynamic damping effects, the aerodynamic damping term is dropped 
from equation (l2) and the aerodynamic moment / M a j- is given by 
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|M a | = [r] l^measured}* Equation (.12) is then applied by use of the pre- 
determined matrices to obtain the structural momenta 

If the hinged-blade equations are used with measured airloads where 
co = Q and the equations involve no aerodynamic damping, the last row 
and last column automatically drop out of equations (.11) and (12) when 
the matrix multiplications are performed, and p h is no longer involved 

in the solution. This behavior of the equations is consistent with the 
physical concept of blade rotation in the plane of no flapping or rigid 
mode resonance. 

Moment of total aerodynamic loading as determined from measured 
structural moment .- A requirement for determining the total aerodynamic 
moment is that the blade flapping angle at the root of the flexible 
blade p^ be measured as well as the structural moment. The procedure, 

then, is to apply equation (ll) and drop the aerodynamic damping term 
from the equation because the structural moments are a result of the 
total loads applied to the blade. 

The aerodynamic moment obtained is thus the total aerodynamic 
moment and includes the moment due to all of the aerodynamic damping 
loads . 

Aerodynamic damping loads .- In some of the foregoing applications 
it may be desired to determine the aerodynamic damping load or moment 
which is proportional to the blade total deflection velocity z in 
order to remove these effects from the measured data. The load and 
moment linearly related to the deflection velocity are given by 


{ L da} = *dQ.M [M j { r - h }] (-"4 

[ Ph J 



culations for the lower harmonics, but p^ would probably have to be 

calculated for the higher harmonics because of its extremely small 
magnitude . 
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Determination of blade natural modes and frequencies .- If the blade 
rotating natural modes and frequencies are desired, the procedure for 
their determination is as follows: 

(1) Solve equation ( 15 ) for unit loads applied at the center of 
each of the ten blade elements. 

( 2 ) Use equation (10.) and the results of step (l) and determine the 
deflections due to the unit loads. 

( 3 ) Set up the numbers from step (2) In matrix form, postmultiply 
by the diagonal mass matrix firQ, and follow standard iteration proce- 
dures. (See ref. 6 .) 

The ten-element breakdown of the blade allows ten degrees of free- 
dom and thus permits the calculation of ten natural modes and frequen- 
cies; however, the use of the standard iteration procedure for obtaining 
accurate results above the fourth natural frequency is rather difficult 
because the iterations converge slowly at these higher frequencies. 


Application to Cantilever and Teetering Blades 

The foregoing procedures were based on the hinged blade. When the 
problem is for a cantilever blade, the only change is to use the corre- 
sponding cantilever -blade equation and matrices. The cantilever-blade 
equivalents of the hinged -blade equations (lO) and ( 25 ) are not given 
bub are obtained by making p^ and p equal to zero in these equations 

The cantilever -blade matrices are given in tables I through V as the 
elements above and to the left of the dashed partition lines. 

In an analysis of a teetering blade, the even harmonic loadings are 
considered as applied to a cantilever blade; the odd harmonic loadings 
are considered as applied to the hinged blade. The reasons for the dif- 
ferent considerations for different harmonics is that even harmonic 
loadings are symmetrical loadings, that is, of equal sign on both sides 
of the hub; the odd harmonic loadings are unsymmetric. 


NUMERICAL EXAMPLES 


In order to illustrate the method, it is applied to a cantilever 
and to a zero-offset hinged blade. Both blades were of the same weight 
and stiffness, with the exception that the stiffness of the hinged blade 
was zero at the hinge. The weight and stiffness distributions are given 
in figure 4 . 
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The aerodynamic loading applied to the hinged blade in forward 
flight is shown in figure 5(a) . The loadings shown are plots against 
spanwise station of assumed Z-^ q and Zj_ g coefficients for the first 

harmonic of measured blade airloads as defined by 

= Zj_ Q COS g Bill 

By using the loading of figure 5(a) and the weight and stiffness of the 
blade shown in figure 4, the structural moment coefficient shown in 
figure 5(b) was calculated by means of equation (.12) . The aerodynamic 
damping term of equation (12) was dropped because measured loads already 
included aerodynamic damping effects. The effects of structural damping 
were not included, although the inclusion of structural damping in the 
equations should improve the accuracy of the calculations for the case 
where the blade is near resonance. However, near blade resonance, the 
effects of small errors in the calculated blade weight and stiffness on 
the calculated results are likely to be large, and the problem requires 
extreme accuracy in all components. 

The structural moments calculated by the modal method of reference 5 
for the loading of figure 5(a) a*"® also shown in figure 5(h) • For these 
calculations the first four natural rotating modes and frequencies were 
used in equations of reference 5 • The aerodynamic damping terms were 
dropped from the equations. The mode shapes and natural frequencies 
were obtained by the method previously outlined in the calculation pro- 
cedures . As can be seen in figure 5(h) , the agreement of the results 
obtained by the two methods is very good. 


The aerodynamic loading applied to the cantilever blade was the 
forward -flight second-harmonic loading shown in figure 6(a). The 
loadings shown are plots against spanwise station of assumed q and 

Zg g coefficients for the second harmonic of measured airloads as 
defined by 


Zg — Z^ q Cos 2)j t *+■ Z^ g sin 2ty 

This loading and the weight and stiffness distribution of figure 4 were 
used in equation (17) to calculate the structural moments shown in 
figure 6(b) . Also shown in figure 6(b) are the results obtained by the 
method of reference 3 by using the same aerodynamic load (from fig. 6(a)) 
and four rotating symmetrical modes and natural frequencies. The results 
calculated by the present method and by the method of reference 3 are 
seen in the figure to be in good agreement. 

The structural moments on the rotating hinged and cantilever blades 
caused by 100-pound loads applied at various spanwise stations are given 
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in figures 7 and 8. The moments were obtained by the use of equations (15) 
and ( 18 ). By using the curves for the moment on the cantilever blade 
caused by 100-pound loads at various stations and the aerodynamic loads 
measured at various stations in the hovering test of reference 7> the 
structural moments in hovering were calculated. The calculated moment 
and the blade deflections for the hovering condition are shown in figure 9- 

The various matrices used in the example calculations are shown in 
tables VII to X. The elements above and to the left of the partition 
lines pertain to the cantilever blade. The matrix of all of the given 
elements pertains to the hinged-blade case. 


CONCLUDING REMARKS 


A matrix method has been derived for determining the structural 
moment and deflections of hinged, teetering, and cantilever blades. The 
method avoids any preliminary mode shape, natural frequency, or quasi- 
static calculations and, as a result, is comparatively short and involves 
only standard matrix procedures. 

The equations are well suited to the determination of the blade 
bending moment at various rotor speeds and for numerous loadings. The 
method puts no restriction on the blade deflection shape and thus permits 
the blade to assume more complex shapes that could be simulated by the 
superposition of four natural modes, as is usually employed in a modal 
method of solution of the differential equation for blade bending moment. 
Thus, the present method should give more accurate results than the modal 
solution for the higher harmonics, where the blade deflection and load 
distribution are of a more complex shape. 

The method includes the effects of the primary inertia and damping 
loads on the blade and also incorporates mainly the following assumptions* 
(l) small angle consideration, ( 2 ) step diagram of the blade mass and 
stiffness, ( 5 ) step integration procedures, and (4) steady average aero- 
dynamic considerations in the damping terms. In addition, the method 
does not include the effects of: (l) radial change of mass position 

with blade deflection, ( 2 ) variation of rotor speed with azimuth angle, 
and ( 3 ) any torsional considerations. Some of the effects of the step 
distributions may be removed by increasing the number of blade elements j 
however, for most blades ten elements should be sufficient. 

The inclusion of structural damping in the equations should improve 
the accuracy of the calculations for the case where the blade is near 
resonance. However, near blade resonance, the effects of small errors 
in the calculated blade weight and stiffness on the calculated results 
are likely to be large, and the problem requires extreme accuracy in all 
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components. These unfavorable factors, in some cases, could be somewhat 
offset by experimental determination of the elements of the influence 
coefficient matrix and the mass of the blade elements and by some 

experimental estimate of the structural damping. In cases where the 
blade is not near resonance with one of the applied load frequencies, 
the structural damping term is relatively unimportant and may be dropped 
from the equations for blade bending moment. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , May 7, 1958. 
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APPENDIX 

BASIC INTEGRATING MATRICES 


Moment of an Arbitrary Loading 
The general bending-moment equation (eq. (l)) is given as 



r R 

pR 



/ l V[ r - r n) dr " 

/ w* 

- z n ) dr 


J r n 7 

J x n y 

/ 


In this section the matrix equivalent of the expression 


M rn - 



(Al) 


is derived. 

The integration of equation (Al) may be accomplished by using a 
procedure based on step, trapezoidal, or parabolic representation of the 
loading. In the present method the step diagram is used, and the load 
is considered to act at the center of each step. 

The procedure for obtaining the moment of a loading is as follows: 

(l) The loading is replaced by a system of equivalent concentrated 
loads as shown in the following sketch: 
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( 2 ) The moment due to the loading is desired at the stations shown 
in the following sketch: 



Arbitrary moment station arrangement 


The loads are centered for 10 equal stations at 0.05(R - h), 0.15(R - h ) , 
0.25(R - h) , and so forth, whereas the moments are taken at stations 
0.85(R - h), 0.75(R - h), and so forth except near the root where stations 
at 0.0125(R - h) and 0.0625(R - h) are used to include possible hub 
effects . 

With the arrangement just shown the moment equation is given as 



In equation (A2) any elements of negative sign resulting from the 
subtraction are to be replaced by zeroes since the moment contribution 
of the loads inboard of a moment station is zero. 
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If the rectangular matrix of equation (A2) is symbolized by [r], 
then equation (A2) may be written as 

{m} = [r] {l| (A3) 

which is the matrix equivalent of 


*v( r - r n)^ 

where the integral is to be evaluated at a selected number of stations. 

The [r] matrix is given in table V for the ten load and moment 
stations of figures 2 and 3- The matrix of table V may be extended to 
include more stations by expanding the equation |m| = [r] {l} as shown 
in equation (A4) , where the rectangular matrix postmultiplied by the 
diagonal matrix is equal to the [r ] matrix-. 
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To extend the matrix to more stations, continue the given number 
sequence of the rectangular matrix and obtain more rows at the top and 
more columns on the left. Keep unity numbers on the diagonal and zeroes 
on the right. The moment stations are then the given stations (shown in 
table V) multiplied by 10/ne, (where n e is the new number of blade 

elements), and the obtained sequence is continued columnwise in order to 
provide additional stations. For example, the elements of the first 
column of a 20-station arrangement would be 1, 2, 3, 4, 5, 6, . . . l6, 
17, 18.87, 19-37, and 19-5(r - h)/(R - h); the load stations would be 
0.025, O.O75, 0.125, O.175, • • • O.925, and 0.975(r - h)/(R - h); and 
the moment stations would be 0, 0.00625, 0.03125, O.O725, 0.125, 0.175, 

. . . 0.875, and 0.925(r - h)/(R - h) . ' 


Moment Due to Horizontal Forces Acting on a Deflected Blade 


The general bending -moment equation of the text (eq. (l)) is 



r R , 


r R / 

M rn = 

X M* 

x n 

• - r n)dr - 

J l E{ z - z n)<^ 

r n 


In this section the matrix equivalent of the expression 


Mm = - %(z - z n )dr (A5) 


is derived. 

The system of equivalent horizontal concentrated loads, as shown in 
the following sketch, is again used, although some other representation 
might produce slightly more accurate results at the expense of simplicity. 
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At station the moment due to horizontal centrifugal forces may- 

be -written as 

M rl =- L H4 ( z h- _ z i) + l H 3( z 5 “ z l) + l H 2( z 2 ~ z l) CA6') 

Collecting terms and -writing equation (A6) in matrix form yields 



n=4 

- Lh4 Lh3 Lh2 - XT L H 



(AT) 


If this procedure is now applied and extended to obtain the moments due 
to the centrifugal inertia forces at ten selected moment stations on the 
blade, the following matrix equation results: 
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where 


*10 = m 10 r 10 = 


L H10 

n 2 


F g = m 9 r 9 



F 1 = m l r l 


%! 



If the rectangular matrix of equation (A8) is symbolized by [f] 

(given as table II), then equation (A8) may be written in shorter form 
as 



(A9) 


r R 

which is the equation M rn = -J Zg^z “ Znjdr in matrix form and is 

r n 

to be evaluated at selected stations. The moments M -m are the moments 

on the blade caused by the centrifugal forces acting on the deflected 
blade . 


In order to extend the matrix equation (A8) to n stations, make 
all of the negative s umma tions on the shown diagonal go from F n to Fj^, 

where k is the subscript of the F values in the preceding column as 
illustrated. The elements off the diagonal of summations are obtained 
by extending the indicated sequence to F n - The new moment and deflec- 
tion stations are 10/n e times the given stations, and the obtained 

sequence is extended columnwise for additional stations. The new blade- 
element arrangement has to have equal-length increments £r . 
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Blade Deflection Due to Structural Moments 

For the blade the general bending-moment equation of the text 
(eq. (7)) is given as 


[l + ig] {M s } = [r] {L a } + ifi co [r] f A] {z} + <o 2 [r] [>] {z} - ^[f] {z} 


In subsequent derivations, a matrix expression for the deflection |z j 

is required. The deflection jzj is composed of a component due to 

the flexible -blade root slope and a component due to blade bending. In 
this section a matrix expression for the deflection component due to 
blade bending is derived. The derivation involves establishing the 
matrix equivalent of the analytical expression 


z r = f f i- dr dr (AlO) 

JO Jo El 

The derivation of the required matrix expression is Illustrated herewith. 
The integration is first illustrated with four stations and then extended 
to more stations. 

The blade stiffness and structural moment distributions are repre- 
sented by the following step diagrams: 



Blade stiffness distribution 
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L_ 

^ .1 

Ar 

Ar 

Ar 


(R - h) I 

(R - h) 

‘ (R-h) ‘ 

* (R - h) ' 


Mi 

1 m 2 

M 5 






M 4 












M 5 







U 1 


.025 .1 .2 .3 


(r - h)/(R - h) 

Blade moment distribution 


The segmental arrangement just noted could be made differently, but 
the particular arrangement shown was chosen for the reasons that; (l) it 
permits a large change of blade stiffness at the blade hub, and (2) it 
gives a moment station where there are large changes in moments when the 
methods are applied to the hinged blade. 


The blade angles 
giyen by the equation 


p r resulting from the bending deflections are 
*r __ / r 


Pr 


= P M_ 

J 0 El 


dr 


[ or the eq. p r = > -M- Ar 

L 0 si 


and can be deduced, for example, from the foregoing distributions as 
follows : 


At 


(r - h) 

(R - h) 


0.025, 


P r = 


M 1 /at\ 
EI-^ / 


At 


(r - h) 
(R - h) 


0 . 10 , 




+ 


m 2 

iij 
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At 


(r - h) 
(R - h) 


0 . 20 , 


Pr = 



At 


(r - h) 
(R - h) 


0.30, 




These distributions may be written in matrix form as follows: 


P .025 


At 

4ei x 

0 

0 

0 


r 

Mi 

P.10 

< 

) = 

Ar 

3 Ar 
4EI 2 

0 

0 


m 2 ' 
> 

P.20 


Ar 

4EI x 

3 Ar 
4EI 2 

Ar 

ei 5 

0 


m 3 

P.30 


Ar 

4EI x 

3 Ar 
4 ei 2 

Ar 

ei 5 

Ar 

EI 4 


M 4 

^ 1 


(All) 
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In order to extend the procedure to obtain the deflections in 
terms of the moments, the following sketch proves helpful: 



At 


(r - h) 

(R - h) 


0.05, 


z = 



M 1 

Ell V 2 


At l r -..a , Q.lg, 
(R - h) 


z = 


z 0.05 + Po.io & = 



'Mt Ar 3Mp £r 

+ — £r 


i 4EIn 


teic 


At 


(r - h) 
(R - h) 


0.25, 


z 


z 0.15 + Po.20 ^ - 
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Extending and collecting terms and -writing in matrix form yields the 
following: 


f 

Lf\ 

O 

ISJ 

V 


—2 

dr 

8EI X 

z -15 


-r yr™2 

3 dr 


£ 

H 

H 

Lf\ 

CVJ 

N 

“V 

> — 

5 dr 2 


8EI! 

z.35 


—2 
7 dr 




3 dr 
4EI; 


6 dr* 


9 dr 
4EI 0 


Symbolizing the square matrix by 




— 



0 

0 


Ml 

0 

0 


m 2 

— 2 
dr 

EIj 

0 

<1 

M3 

— 2 

2 dr 

EI3 

dr 2 

EI4 


M4. 

M - 


gives equation (A12) as 


(A12) 




(A13) 


Rearranging equation (A12) to have the blade tip values in the 
upper rows and extending to obtain deflection values at the twelve 
selected stations results in the following matrix. (See figs. 1, 2, 
and 3 for station location and the definition of the elements.) 
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(A14) 



52 


MCA TN 4304 


Equation (Al4) may be written as follows: 
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To extend to more stations , simply extend the established number 
sequence in the rectangular matrix. Keep unit values on the extended 
diagonal and zeroes below. The procedure for obtaining the new moment 
and deflection stations is as previously discussed in the jjr J and 
JV] matrix derivations. The values of El must pertain to the new 
moment stations. 

The results obtained by using either equation (Al4) or (AI 5 ) have 
been compared with results obtained by graphical numerical double 
integration of a known M and El distribution, and the results 
agree very well. If desired, the elements of equation (Al4) or (AI 5 ) 
could be determined experimentally from the actual blade. 
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TABLE I,- MODIFIED IDENTITY MATRIX 


fr-°aj{o} 

LoJ | o 





Moment 

station^ 

(r 

- h)/(R 

- h) 



i 

J 0 


Moment 
station, 
r - h 
R ~ h 

0.85 

0.75 

0.65 

0.55 

0.45 

0.35 

0.25 

0.15 

0.0625 0.0125 

, p h 

1 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 0 
1 


0.85 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 0 
1 


.75 

0 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

1 

1 0 
1 


.65 

0 

0 

0 

1.0 

0 

0 

0 

0 

0 

0 

1 

1 0 
1 


.55 

0 

0 

0 

0 

1.0 

0 

0 

0 

0 

0 

1 0 

i 


.45 

0 

0 

0 

0 

0 

1.0 

0 

0 

0 

0 

I 

1 0 
1 


.35 

0 

0 

0 

0 

0 

0 

1.0 

0 

0 

0 

1 

1 0 
1 


.25 

0 

0 

0 

0 

0 

0 

0 

1.0 

0 

0 

1 

i 0 
l 


.15 

' 0 

0 

0 

0 

0 

0 

0 

0 

1.0 

0 

1 

1 0 
1 


.0625 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1.0 

1 

* 0 
1 


.0125 

— 

- 

— 

— 

— 

— 

— 

— 

— 


-J — 


— 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 0 


0 
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TABLE II.- CENTRIFUGAL-FORCE MATRIX [f] 


[ F n = Vn] 






Deflection station 

, (r 

- h)/(B - 






Moment 

0.95 

0.85 

0.75 

0.65 

0.55 

0A5 

0.35 

O.25 

0.15 

O.O625 

0.05 

O.OI25 


station, 
r. - h 
R - h 


- F 10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 ^ 


0.85 

F xo 

f 9 

n=10 

- 5 Z 

n=9 

F n 0 

0 

0 

0 

0 

0 

0 

0 

0 


•75 

F 10 

f 9 

f 8 

n=10 

n-8 

F n 0 

0 

0 

0 

0 

0 

0 

0 


.65 

F 10 

F 9 

f 8 

F 7 

n=10 

-XI F n 

n =7 

0 

0 

0 

0 

0 

0 

0 


.55 

F 10 

f 9 

f 8 

f 7 

f 6 

n=10 

"S Fn 

0 

0 

0 

0 

0 

0 


.45 

F 10 

F 9 

f 8 

F 7 

f 6 

F 5 

n=10 

-2= 

n =5 

F n 0 

0 

0 

0 

0 


•35 

F 10' 

F 9 

f 8 

F 7 

f 6 

F 5 

f 4 

n=10 

n =4 

L 0 

0 

0 

0 


•25 

F 10 

F 9 

f 8 

F 7 

f 6 

F 5 

f 4 

F 5 

n=10 

n =3 

F n 0 

0 

0 


• 15 

F xo 

F 9 

f 8 

F 7 

f 6 

F 5 

*k 

F 3 

F 2 

it=10 

"XT F n 

n=2 

0 

0 


.0625 

F 10 

F 9 

f 8 

F 7 

f 6 

F 5 

f 4 

F 3 

F 2 

0 

F 1 

n=10 

-ST 

n=l 


.0125 

_ F 10 

F 9 

V 

F 7 

f 6 

F 5 

F 4 

F 3 

F 2 

0 

y i 

0 


0 
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TABLE III.- DEFLECTION MATRIX 


.[%K r - *0 






Moment 

station. 

(r - h)/(R 

- h) 


j „ 

1 

0.85 

0.75 

0.65 

0.55 

0.45 

0.35 

0.25 

0.15 

0.0625 

1 % 
0.0125 \ 

1 

EI 10 

2 Sr 2 
Elg 

3 2SF 2 
Elg 

4 5^ 
EI 7 

5 S 2 
ei 6 

6 ^ 
EI5 

ei 4 

8'$* 

EI 5 

27 AF 2 
4EIg 

\ (r " h) o*95 

i 

0 

5 2 

ei 9 

— 2 

2 Ar 

EI B 

5 a? 2 

4 

Elg 

5 g 5 
EI 5 

6 S 2 
EI k 

7 

ei 5 

24 Sc^ 
4EIg 

^ | (r - h) °- 8 5 
1 

0 

0 

z ? 2 

ei 8 

2^ 

El? 

5 A 2 
ei 6 

4S 2 

EI 5 

5 Zr 2 
EI4 

6 Ar 2 
ei 5 

el S 2 

4EI2 

i 

0 

0 

0 

s 2 

EI-jr 

2 A? 2 

5 Ar^ 
El^ 

4S 2 

EI 4 

5 A 2 
EI 3 

18 At 2 
4EI2 

i 

* <*■ - ^0.65 
1 

0 

0 

0 

0 

Sr 2 

ei 6 

2^ 

ei 5 

5 Aiv 

ei 4 

4 ZF 2 
ei 5 

2 • 

15 Zr 

4EI2 

1 

0 

0 

0 

0 

0 

ei 5 

2 ^ 
ei 4 

is: 

EI 5 

12 Ac^ 
4EI2 

1 

< r -«0A 5 

1 

0 

0 

0 

0 

0 

0 

A^ 

EI4 

2 A? 2 
EI 5 

9 ^ 

4EIg 

1 

( r - h io.35 

1 

0 

0 

0 

0 

0 

0 

0 

H — 2 

ei 5 

6 ^ 
4EI 2 

1 

(r “ to) p.25 

1 

0 

0 

0 

0 

0 

0 

0 

0 

5 25^ 
4EI 2 

7-2 1 

1 (r - h)o.i5 
mil l 
l 

0 

0 

0 

0 

0 

0 

0 

0 

5 Zr 2 
32EI 2 

i 

5JSF 2 ) / h i 

h ) 0 .0625 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

* 

| < r - h) o.o5 

l 

0 

0 

0 

0 

0 

0 

0 

0 

0 

^• j (r-h) 0 .oia5 


Deflection 
etatiqn, 
r - h 
R - h 

0.95 

.85 

.75 

.65 

•55 

.45 

*55 

.25 

.15 

.0625 

•05 

.0125 
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TABLE V.- INTEGRATING MATRIX [r] 
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9.5 Ar 8.5 Ar 7.5 Ar 6.5 Ar 5-5 Ar 4.5 Ar 3-5 Ar 2.5 Ar 1.5 Ar 






TABLE VI.- AERODYNAMIC DAMPING MATRIX [V) 
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■ .0125 




NACA TN k 30 k 


kl 


3 



2144.988 1927-752 1700. 957 1113.766 1253.056 789.338 832.589 647.621 1995*864 0 1533*551 




TABLE VIII. - VERTICAL-INERTIA, MATRIX [M V “L lb-sec% FOR EXAMPLE BLADE 
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.503869 .452857 .399562 .261636 .294355 .185423 .195589 .152137 .468867 0 .360292 




TABLE IX. - DEFLECTION MATRIX 2« X -10° i ir - h FOR EXAMPLE BLADE 
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.100 I 1.144 .0125 





TABLE. X.- INTEGRATING MATRIX [r], in., FOR EXAMPLE BLADE 



86.925 77-775 68.625 59-475 50.525 41. 175 32.025 22.875 13-725 0 4.575 
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Figure 1+.- Blade weight and stiffness used in the calculations. 
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Figure 9.- Bending moment and deflections on cantilever blade in hovering. Thrust = 490 lbj 

OR = 497 ft /sec. 



